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Abstract 

We consider the diffuse gamma ray data from Fermi first year observa- 
tions and compare them to the gamma ray fluxes predicted by Dark 
Matter annihilation or decay (both from prompt emission and from In- 
verse Compton Scattering), for different observation regions of the sky 
and a range of Dark Matter masses, annihilation/decay channels and 
Dark Matter galactic profiles. We find that the data exclude large re- 
gions of the Dark Matter parameter space not constrained otherwise 
and discuss possible directions for future improvements. Also, we fur- 
ther constrain Dark Matter interpretations of the e ± PA ME LA/ Fermi 
spectral anomalies, both for the annihilating and the decaying Dark 
Matter case: under very conservative assumptions, only models produc- 
^ ing dominant ly and assuming a cored Dark Matter galactic profile 

can fit the lepton data with masses around ~ 2 TeV. 
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1 Introduction 

A wealth of astrophysical and cosmological data have revealed the crucial gravitational role of 
otherwise essentially non-interacting particles, dubbed Dark Matter (DM). This provides one of the 
most tantalizing hints for physics beyond the Standard Model (SM). However, in order to identify 
the nature of these particles, one needs information on their mass and interaction properties. 
One strategy is to explore the possible indirect signatures coming from their Standard Model 
annihilation (or decay) products in the halo of our Galaxy, or beyond. For typical candidates in 
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the category of weakly interacting massive particles (WIMPs) these products fall in the realm of 
investigation of high-energy astrophysics (typically E > 1 GeV), as gamma-rays, neutrinos, etc. 

The Fermi Gamma-ray Space Telescope [1] team has recently released its first year data. 
While the analyses for their detailed interpretation in terms of known astrophysical sources are 
still ongoing, preliminary spectra from different regions of the sky have been presented at several 
conferences, most notably the Fermi Symposium [2]. These spectra extend to comparable or 
higher energies than the predecessor mission EGRET, but above all have significantly reduced 
systematics. In particular, "anomalous" (with respect to simple propagation models predictions) 
bumps of radiation above the GeV have not been confirmed, see e.g. [3]. 

Given these new data, it is timely to provide a first assessment of their power in constraining 
DM properties. Rather than considering a specific scenario beyond the Standard Model, we adopt 
a phenomenological approach, namely: i) assume DM annihilation or decay in a single dominant 
channel (e + e~, /i + yU~, r + r~, W + W~ , bb or tt); ii) for different choices of the galactic DM halo 
profile, derive constraints in the m x -(av) or m x -T& cc plane, for the cases of self-annihilating 
or decaying Dark Matter respectively. Here m x is the mass of the DM particle \i k av ) its self- 
annihilation cross section and Td ec its half life (for the decay case, we limit to the leptonic channels, 
as discussed below). 

In Sec. 2 we deal with a short description of the formalism, the datasets used and the choices 
for the galactic Dark Matter halo profile, while the constraints that we obtain are reported in 
Sec. 3. In deriving our bounds we assume no prior information from other cosmic ray data, nor 
we worry about astrophysical sources of gamma rays which are surely present and probably even 
dominating the signal. This is clearly an over conservative approach: in Sec. 4 we comment on 
the expected reach in exploring the DM parameter space after a better understanding of the 
astrophysical background is achieved. 

Since, motivated by peculiar spectral features, a plethora of models have appeared trying to fit 
the positron fraction and total electron data from PAMELA, Fermi and HESS [4, 5, 6, 7] in terms 
of Dark Matter (as we review in Sec. 2.3), in Sec. 3 we also comment upon the consistency of these 
scenarios in the light of the new diffuse gamma ray data, as a specific application of the general 
constraints derived above. We limit ourselves to modes of annihilation in a pair of charged lepton 
final states at tree level since, on one hand, they appear the only ones not excluded by other 
constraints (most notably antiproton fraction [8, 9]), on the other hand they are significantly less 
model-dependent than scenarios with annihilations proceeding via additional, non SM states [10]. 
For the case of decaying DM, which has been invoked as alternative explanation of the lepton 
anomalies, a fortiori we limit ourselves to leptonic channels. We anticipate that explanations 
based on DM models are mostly excluded as leading contributors to the PAMELA/ Fermi features. 
The only exceptions (but see discussion in Sec. 4) appear to be for models assuming mostly 
final states and, for annihilating models, also requires cored galactic DM profiles. Note that even 
models producing dominantly r ± , previously allowed, are now excluded (or strongly disfavoured 
for sufficiently cored profiles) by the data: this applies also to decaying scenarios. 

A discussion of the caveats and further constraints, perspectives, etc. is reported in Sec. 4, 
together with our conclusions. 

2 Datasets and constraints 
2.1 Calculation of the signal 

The differential flux of photons from (self-conjugated) dark matter annihilations in the galactic 
halo, observed from a given direction in the sky, is written as 

d$|X 1 M & f 1 f p^o[r(sM Y ^ 

dedQ 2 4 7r Q mlL aS r & \ p & ) de ' {l) 
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where the coordinate r, centered on the Galactic Center, reads r(s, ip) = (r + s 2 — 2 r s cos tp) 1 / 2 , 
r = 8.33 kpc [11] is the most likely distance from the Sun to the Galactic Center (GC) and ip 
is the angle between the direction of observation in the sky and the GC. In terms of the galactic 
latitude b and longitude I, one has cosip = cos b cos / . The coordinate s parameterizes the distance 
from the Sun along the line-of-sight (los). The gamma ray energy is denoted by e. Particle physics 
enters via the DM mass m x , the annihilation cross section (av), and the photon differential energy 
spectrum dN/de per annihilation. 

This spectrum is in general constituted of two main components: i) the prompt gamma 
rays, originating from the fragmentation of the primary products of annihilation (essentially via 
bremsstrahlung of charged particles in the shower and production of tt° and their subsequent de- 
cay into 77), and ii) the Inverse Compton (IC) gamma rays, produced by the upscattering of the 
low energy photons of the starlight, the infrared light and the CMB (denoted collectively as the 
InterStellar Radiation Field, ISRF) by the energetic e injected by Dark Matter annihilations. 1 
We compute the IC gamma ray spectrum as discussed in detail in [12], adopting in particular 
the approximation of neglecting the diffusion processes of the e injected by Dark Matter. While 
this can have an impact on small internal galactic regions (where the DM density has the highest 
gradient), this is instead a very good approximation for large observational regions away from the 
GC (amounting to a correction of less than a factor of 2 on the amplitude of the flux already for 
intermediate latitudes [12, 13]). The total spectrum dN/de is strictly independent of direction 
only for the prompt component, while for the IC it carries a dependence on the characteristics of 
the background light on which the IC scattering occurs. We will however work under the approx- 
imation of constant ISRF over each of the observation regions that we consider (see below), as 
already done in [12]; whenever the region is away from the GC, the dominant background for IC 
processes is in any case the Cosmic Microwave Background (CMB), which is obviously identical 
in each region. The flux from a region AQ of the sky therefore simply rearranges in terms of the 
usual average geometrical factor J ann as 



de 2 4tt m 2 x de J AQ J los r Q \ p Q 

In the case of decaying dark matter one has 



de dQ 4vr & m x J los r \ p / de 



and thus (the same discussions as above apply) 

d$ halo _ r PQ J AO diV T AO_ [ ^O^n [ ds ^Phalo[r(s, 



de 4vr m x de J An J los r Q \ p Q J 

where T = r^ 1 is the decay rate (inverse lifetime) and the spectrum now refers to the photons 
generated in each decay process, including both the prompt and the Inverse Compton emission. 

For the decaying DM scenario, we will also need the cosmological flux, i.e. the flux of gamma 
rays due to cosmological DM decays integrated over redshift. As we will see, it leads to powerful 
constraint and — differently from the analogous annihilating DM signal — it does not depend on 



1 Othcr contributions, such as a line component from loop-suppressed annihilations into 77, jZ or 7/1 (with Z 
and h the Standard Model gauge and Higgs boson) and Internal Bremsstrahlung gamma rays from the charged 
particles in the annihilation diagrams could be present, but, as they are model dependent, we do not consider them 
here. In general they introduce features in the total gamma ray spectrum that would make stronger the constraints 
from an analysis like ours. 
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halo profiles or DM substructures, relying only on robustly known quantities. It is straightforward 
to show that the differential spectrum in Earth-measured energy e from the entire sky writes 



d$ 



dec 



n 



de 



DM Pc,0 



dz 



-r(e(z),z) 

Hiz) lie 1 



e(z),z) 



(5) 



where the Hubble function, H(z), writes explicitly H(z) = Hq^JVLm{)- + z) 3 + £Ia where H is 
the present Hubble expansion rate and Qm and are respectively the matter and cosmological 
constant energy density in units of the critical density, p Ci o- This flux originates from the integral 
along the time of travel of a photon emitted in the past, when converted as an integral in redshift 
z according to the standard relation dt/dz = 1/((1 + z) H(z))\ the (1 + z) factor is canceled by 
an identical (1 + z) at the numerator due to the redshift in the energy of the same photon. In 
turn, the rescaling factor with redshift of the DM density ((1 + z) 3 ) cancels with an identical 
factor coming from the volume expansion of the Universe. Finally, the factor e~ T ^ e ' z ^ accounts 
for the finite optical depth, r, of the universe to high energy gamma rays due to scattering with 
the extragalactic background light (see e.g. [14] for a more detailed discussion). For energies 
e < 100 GeV, the latter has only a moderate impact. The gamma ray spectrum dN/de is again 
the sum of the prompt (P) and IC contributions, as produced at any redshift z. The dependence 
on z in the prompt component is just due to the redshifting of the energy in the argument, i.e. 



diV7de (e(z) 



dN /de(e(l + z),0). For the IC component, the ISRF consists in this case 



of CMB photons only, rescaled in energy according to the redshift. The integral over redshift is 
dominated by the range up to 2 ~ 20. We notice that, in the (well fulfilled) limit where the IC 
loss time is much faster than the Hubble time, the electrons lose all their energy via IC, and all 
scatterings happen in the Thomson regime, one has 



d7V IC 
de 



[eizuz 



1 diV IC 
1 + z de 



e,0). 



(6) 



The above formula can be also understood heuristically as follows: the IC cross section in the 
Thomson regime is energy- independent; since the energy spectrum of the injected particles is the 
same at any z, the IC photon at the production epoch z will have an energy proportional to the one 
of the upscattered background photon. The latter is {1 + z) times the current one, which is exactly 
the factor compensated by the subsequent redshifting. So, the IC spectrum is universal and equal 
to the one calculated at z — 0. The remaining factor (1 + z)~ x can be immediately understood 
since dt = dz/(H(z)(l + z)): apart for the absorption correction, the integral in Eq. (5) for the 
IC is the "spectrum per decay event" times the number density f2 DM p c $/m x of the metastable 
particles, times the probability of decay over the lifetime of the universe tjj, i.e. tjjT. 

The Fermi collaboration has presented preliminary results on the (maximal) residual, isotropic 
gamma-ray flux present in their data (see below). When using these diffuse data to constrain 
decaying DM properties, one has to account for the fact that the DM decays contribute with two 
terms to the above flux: i) a truly cosmological flux, which is of course isotropic at leading order, 
ii) the residual emission at high latitudes from the DM halo of our Galaxy: for the approximation 
of spherical halo, this is minimum in the anti-GC and grows very slowly towards the intermediate 
latitudes (see e.g. [15]). In formulae, the predicted DM flux that we will compare with Fermi 
isotropic diffuse 7-ray data is 



d$ dcc 

isotropic 

de 



d$ 



dec 



de 



+ Air 



dedfi 



(7) 



anti-GC 



For typical DM decay channels and for any DM halo profile, we find that the two contributions are 
of comparable amplitude. Note that a similar equation would hold for the annihilating DM case 
with two significant differences: a) the halo signal has a stronger dependence from the angular 
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distance from the GC; b) while for the decaying case the cosmological signal is readily calculated 
in terms of known quantities, see Eq. (5), for the annihilating case a dependence is introduced on 
the halo profiles and the dumpiness of DM halos: this is the reason why we do not include this 
constraint in the following, but in Sec. 4 we shall discuss its likely importance for cored halos. For 
a comparison of the role of the two terms in the annihilating case, see also [16]. 

For the galactic distribution of Dark Matter we consider the cases of a Navarro, Frenk and 
White (NFW) [17], Einasto [18, 19] and cored Isothermal [20] profiles. The first (peaked as r -1 
at the GC) is a traditional benchmark choice motivated by N-body simulations, the second (not 
converging to a power law at the GC and more chubby than NFW at the location of the Sun) 
is emerging as a better fit to more recent numerical simulations. Cored profiles, such as the 
Isothermal one or the Burkert profile [21], might be instead more motivated by the observations of 
galactic rotation curves, but seem to run into conflict with the results of numerical simulations. As 
long as a convergent determination of the actual DM profile is not reached, it is worth considering 
a range of possible choices. The functional forms of these profiles read 

T f T \ ^ 

PNFwO") = Ps— 1 + — , (8) 

r \ r s / 



pEm{r) = p s exp 



0.17 



0.17 

- 1 



(9) 



Piso T (r) = P \ 2 - (10) 
1 + {r/r s ) 

A separate but related issue is the one of fixing the parameters r s and p s that enter in each 
of these forms. We apply the following criterion. We fix the NFW parameters at the benchmark 
values r s = 20 kpc and p s = 0.26. This produces a DM halo with a density of p = 0.3 GeV/cm 3 
and with a total DM mass contained in 50 kpc (i.e. within the distance of the Large Magellan 
Cloud) of M 50 = 3.6 x 1O U M . We then impose that the other profiles obey the same constraints 
on p Q and M50. This is motivated by the following considerations: on one hand, the total mass 
estimate within 50 kpc is quite robust [22], on the other hand recent analyses suggest that the 
best-fit value for the density of the DM at the solar distance is the same for different models to 
within ~ 7% [23]. The parameters that we adopt are thus given explicitly by: 



DM halo model 


r s in kpc 


p s in GeV/cm 3 


NFW 


20 


0.26 


Einasto 


21.8 


0.05 


Isothermal 


3.2 


2.31 



We emphasize that in no way we are claiming that the choices for the parameters {p s ,M 50 } are 
the optimal ones, but by the above criteria we are "normalizing" the alternative profiles in a more 
physically motivated way, commensurate with the benchmark profile which is often employed in 
the literature. In any case, we note that employing slightly different normalizations (such as those 
already employed in [12]) has an overall very small impact on the analysis. 



2.2 Fermi diffuse gamma-ray observations 

The Fermi collaboration has presented results on diffuse gamma rays from the first year analysis 
in terms of data points in several different observational windows. For the annihilating Dark 
Matter case, we focus in particular on inner galaxy and intermediate latitude regions. Namely, 
we consider: 
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Region 


latitude b & 
longitude I 


IsoT 


J arm 

NFW 


Einasto 


IsoT 


"dec 

NFW 


Einasto 


ISRF 


'3x3' 


0° < \b\ < 3° 
0° < l < 3° 


35 


374 


635 








GC 


'5x30' 


0° < |6| < 5° 
330° < I < 360° 
0° < I < 30° 


21 


5G 


93 








GC 


'10-20' 


10° < \b\ < 20° 
0° < I < 360° 


3.35 


3.46 


4.23 


2.45 


2.38 


2.47 


h = 5 kpc 


'Gal Poles' 


60° < |6| < 90° 
0° < I < 360° 


0.92 


0.96 


0.94 


1.73 


1.69 


1.67 


CMB 



Table 1: Summary of the observational regions that we consider, with the corresponding values of 
the average J factor for different DM halo profiles and the adopted Inter-Stellar Radiation Field. 

o A square region that includes the Galactic Center, at galactic latitude 0° < |6| < 3° and 
galactic longitude 357° < / < 360° and 0° < / < 3°, that we denote as '3° x 3°' region. We 
use the preliminary Fermi data points as presented in [24]. 

O A rectangular region that also includes the Galactic Center, at galactic latitude 0° < \b\ < 5° 
and galactic longitude 330° < I < 360° and 0° < I < 30°, denoted as '5° x 30°' region [24]. 

O The intermediate latitude region defined by 10° < |fe| < 20° in galactic latitude, at all 
galactic longitudes (denoted in the following as '10° — 20° strips'). The released Fermi data 
points (from [3]) improve on the former release in [25, 26], confirming the results. 

o The high latitude 'Galactic Poles' region, identified by |fe| > 60° at all longitudes. We use 
the data points from [27]. 

For the decaying Dark Matter case we consider the two intermediate/high latitude galactic 
regions discussed above ('10° — 20° strips' and 'Galactic Poles' ) and, in addition: 

o The measurements of the isotropic diffuse flux presented in [28] (see [29] for a similar kind 
of analysis, limited to decaying DM models). 

In all the regions, the preliminary Fermi data points that we use extend to about 100 GeV. 
Ongoing work in the collaboration will likely extend this range to higher energies. We comment on 
the impact of this below. The uncertainties that we use are generally quoted to include statistical 
and systematic errors, but a more detailed analysis is in the making by the collaboration. 

In order to compute the ICS signal from DM annihilations or decays, we model the ISRF 
for each region as discussed in [12], in terms of renormalized black body spectra that reproduce 
the background light characteristic of the different areas. For the inner galaxy regions ('3° x 3°' 
and '5° x 30°') we use the ISRF computed for the Galactic Center (R = 0, h = 0, in galactic 
cylindrical coordinates) in [30], that features an important contribution at small wavelengths of 
starlight photons, because of the higher density of stars in the galactic bulge. For the '10° — 20° 
strips' we use the ISRF computed for distances above the galactic plane, h — 5 kpc, in [31]. For 
the Galactic Poles and of course for the cosmological flux we use the CMB only. 

2.3 Details on the DM interpretations of the e ± anomalies 

As mentioned in the Introduction, the anomalous PAMELA, Fermi and HESS data have been 
interpreted in terms of DM annihilations [8, 13] or decay [32]. We recall here briefly the main 
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features of the experimental data and of their DM interpretations, without entering in the details 
of any specific particle physics model. 

The PAMELA satellite has reported a rise of the positron fraction e + /(e + + e~) above the 
expected declining astrophysical background from ~10 GeV up to at least 100 GeV [4]. At the 
same time, the ratio p/p fluxes is consistent with expectations based on astrophysical secondaries, 
up to the maximal probed energy of about 200 GeV [34, 35]. Fermi [5, 33] has reported a e + + e~ 
spectrum suggestive of an additional harder component on the top of a smooth astrophysical 
spectrum, the latter falling slightly more steeply than E~ 3 . The HESS Cerenkov telescope, too, 
has published data [6, 7] in the range of energy from 600 GeV up to a few TeV, showing a power 
law spectrum in agreement with the one from Fermi and eventually a steepening at energies of a 
few TeV. 

We perform the fits to these data 2 with the use of DM generated e + , e + + e~ and p spectra, as 
discussed in detail in [8] : we find the best fit values by scanning over the propagation parameters of 
charged cosmic rays and over the uncertainties on the slope and normalization of the astrophysical 
electron, positron and antiproton background. Our analysis is again conservative in the fact 
that we consider the slopes and normalizations of the backgrounds of the different species as 
uncorrelated, while a significant correlation exists (e.g. both positrons and p are produced in the 
same collisions as astrophysical secondaries). A proper treatment of this effect would reduce the 
"allowed parameter space" for DM models. We do not include any galactic boost factor due to 
substructures within the Milky Way halo. We also stress again that we consider only the case in 
which the annihilation or decay proceeds directly into a SM particle-antiparticle pair. I.e. we do 
not address here the cases of one-step annihilation models [10], of decay into three leptons (see 
e.g. [39]) etc. 

The resulting allowed regions on the plane m x —(av) are shown in fig.s 2 and 3 for the annihi- 
lating DM case and on the plane m x -T^ ec in fig. 4 for the decaying DM one. The regions that allow 
to fit the PAMELA data alone (for positrons and antiprotons) are individuated by green and yellow 
bands, for 95% C.L. and 99.999% C.L., corresponding to Ax 2 {m x , (crv)) = y^iyn^ (av))— xLw — 6 
and Ax 2 — 23, with 2 d.o.f. (analogously for decay). The addition of the Fermi and HESS data 
allows to individuate a narrow range for the mass of the DM particle, around a few TeV. It also 
selects a DM that annihilates or decays into leptons only, since, for such a mass, the excess in 
antiprotons that would arise from non-leptonic channels is excluded by PAMELA p data. The 
regions that fit PAMELA+Fermi+HESS combined are represented in the figures by red and orange 
areas (for the same confidence levels as above). Notice that the smoothness of the Fermi spectrum 
forbids a reasonable fit with the DM DM — > e + e~ or DM — > e + e~ channel, that would produce 
too peaked a feature: the reduced (Xmm)red turns out to be well above 2 for all DM profiles, so 
we plot no allowed region. The typical annihilation cross sections that are required are of or- 
der 10~ 23 cm 3 /sec up to 10 _20 cm 3 /sec or more, depending on the mass of the candidate and the 
annihilation channel. For decay, t& cc w 10 26 sec is needed. 

2 We consider the PAMELA data for positrons at energies larger than 10 GeV only, where the uncertainty due 
to solar modulation is not present. We use the new preliminary analysis of p data from PAMELA, as presented 
e.g. in [35], that extends to about ~ 200 GeV and provides somewhat reduced uncertainty bars. We include the 
systematic error band on the Fermi and HESS datapoints. We add the low energy datapoints in the e + + e~ from 
Fermi as presented in [33]: this however does not change the fit regions much. We do not include in the fit the 
data from the ATIC-2 and 4 [36, 37] and PPB-BETS [38] balloon experiments, that in particular had found a peak 
at about 500-800 GeV (possibly indicative of a DM mass around the TeV scale), as this feature has later been 
questioned by the Fermi findings. Since however Fermi+HESS also pin down the mass in a similar range, our results 
would easily apply also to the case in which ATIC is reincluded. 
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Figure 1: Gamma ray fluxes for a few sample DM candidates, compared to the Fermi datapoints 
in the different observation regions that we consider. See text for details. 



3 Results 

We start by reporting, in Fig. 1, the gamma-ray fluxes for a few typical TeV DM candidates with 
large annihilation cross section in leptonic channels (of the type invoked to explain the anomalies 
in e data), in different angular windows and for different halo profiles. For each case we plot the 
total gamma-ray flux and its different components: the prompt gamma-ray emission and the ICS 
emission on StarLight (SL), on InfraRed light (IR) and on the CMB. Similar plots can be drawn 
for the decaying DM case. 

As apparent, in all these cases the spectral shapes of the curves of the DM signals are very 
different from what is observed. In a e 2 d$/de plot, the Fermi data point show a decreasing 
behaviour (that often seems well-accounted for by a simple, likely astrophysical, power-law), while 
a curve rising up to ~ TeV energy, possibly with a "double bump" feature (characteristic of the high 
energy prompt and low energy ICS emissions) is expected from DM. This immediately reasserts 
that a significant astrophysical signal is needed to account for the data, confirming the conservative 
approach of our analysis. 

The first two panels of Fig. 1 show the predicted signal in the '3° x 3°' and '5° x 30°' regions from 
a DM candidate of mass 1.5 TeV, annihilating with 100% B.R. into fi + n~ with a cross section 
of 3 10 -23 cm 3 /sec, assuming an NFW or Einasto (i.e. those suggested by numerical N-body 
simulations) respectively. It is evident that the predicted signal overshoots the data points, very 
evidently in the first case but also significantly in the second case. These kind of DM candidates 
are therefore clearly excluded by observations. 
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The third panel of Fig. 1 shows the predicted signal in the 'Galactic Poles' region for the same 
DM candidate, assuming an Isothermal profile. In this case the signal lies well below the Fermi 
data points. This kind of scenario cannot be excluded or explored yet in this observational region. 

The second row of panels in Fig. 1 shows the current situation for the Isothermal cored profile 
choosing a 3 TeV DM annihilating into t + t~ with a cross section of 2 10~ 22 cm 3 /sec. The Fermi 
data points in each region are in clear tension with the predicted DM signal. Note that for the 
purposes of the following plots we shall not depict these models as "firmly excluded" , although this 
would clearly be the conclusion of a proper likelihood analysis combining the different channels 
(not to speak of further constraints from non-7 channels). 

As an aside, the plots in Fig. 1 also allow to appreciate the different spectral features of the 
predicted DM signals originating from the different ISRF in the different regions. For the inner 
regions ('3° x 3°' and '5° x 30°') the ICS signal is predicted to be dominated by the contribution 
on StarLight. For the intermediate region ('10° — 20° strips') the lower energy contribution on 
CMB is more important. At the 'Galactic Poles' the CMB contribution only is considered: the 
other components are subdominant. The prompt contribution is particularly relevant for the t + t~ 
channel, as it originates from the 7r° — > 77 from the hadronic r decays. At the higher latitudes 
it is well distinct from the IC emission, as the latter is located at lower energies. 

The above examples guide us in the interpretation of the exclusion plots, discussed below. 
We derive the constraints on the annihilation cross section or the decay half-life by the following 
conservative prescription. For each observational region with the corresponding data points, fixed 
an annihilation channel, a DM distribution profile and a DM mass, we impose that the ICS 
signal must not exceed any of the experimental data points by more than 3a. This determines 
a maximum annihilation cross section or a minimum half life. Notice once again that we do not 
assume anything on the astrophysical gamma ray background, although it must be clearly present 
and dominant by visual inspection of the plots. We will return on this point in Sec. 4. 
In Fig. 2, we report the exclusion plots for leptonic annihilation modes. We find that: 

o For profiles of the NFW or Einasto type, i.e. those suggested by numerical N-body sim- 
ulations, current data from the inner Galaxy exclude not only DM scenarios explaining 
simultaneously Fermi+HESS and PAMELA features, but also PAMELA data alone to a high 
confidence level. The only small region left in parameter space (and with a low confidence 
level) is the one at m x just above 100 GeV, i.e. just above the highest energy point pre- 
sented by PAMELA. This implies a final test at AMS-02 of the DM origin of the PAMELA 
anomalous positron fraction: should the "excess" continue to rise well above 100 GeV, the 
corresponding DM mass would enter the region already excluded, in clear conflict with the 
Fermi diffuse gamma ray data analyzed here. 

o For "cored" profiles, the data fail to fully exclude, in this simplified analysis, the regions 
explaining PAMELA+Fermi+HESS features. The r channel is however already in strong 
tension with existing data (and likely the 2 3cx best fit region would be excluded to signif- 
icant confidence level by a combined analysis of the 7-channels), while the yir annihilation 
is still allowed by a factor ~ 2 4- 3. 

o For NFW or Einasto profiles, already this extremely conservative analysis gives constraints 
on (crv) which, for typical masses m x ~ 100 GeV, are above the value expected for a S-wave 
annihilating thermal relic by a factor of order 10 or less, especially for the e mode. 

In Fig. 3, we report the exclusion plots for hadronic annihilation modes. We can see that: 

o In the ranges where the different channels are open, the bounds are almost independent from 
the channel: this follows from the quasi-universality of the gamma spectrum arising from 
the fragmentation of sufficiently "heavy" SM particles. 
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Figure 2: T/ie regions on the parameter space m x -(o~v) that are excluded by the diffuse galactic 
gamma ray measurements by the Fermi satellite. The first column of panels refers to DM annihila- 
tions into e + e~ , the second into and the third into t + t~ ; the three rows assume respectively 
an NFW, an Einasto and a cored Isothermal profile. Each panel shows the exclusion contour due 
to Fermi observations of the '3° x 3°' region (blue short dashed line), '5° x 30°' region (orange 
dashed line), the 10° — 20° strip' (red long dashed line) and the 'Galactic Poles' \b\ > 60° region 
(black long dashed line). We also report the regions that allow to fit the PAMELA positron data 
(green and yellow bands, 95 % and 99.999 % C.L. regions) and the PAMELA positron + Fermi and 
HESS data (red and orange blobs, 95% and 99.999% C.L. regions). 



10 



DM DM -» bb, Einasto profile 



DM DM -» WW, Einasto profile 



DM DM -» tt, Einasto profile 



10 _2o p 




10- 



- FERMI 3x3 
FERMI5"x30" 
FERMI 10"- 20" 
FERMI Gal. Poles 

to 3 



to- 24 

1(T J5 

io- m 



m x [GeV] 
DM DM -> bb, NFW profile 




FERMI3'x3" 1 10 " 

FERMI 5' x 30" 

FERMI 10"- 20" | 10- 25 

FERMI Gal. Poles : 

3 io- 2 " 



10 2 I0 3 
m x [GeV] 

DM DM -» bb, Iso profile 



id 4 




10 _20 P 



iir 21 



FERMI 3x3 
FERMI 5" x 30" 
FERMI 10"-20" 
FERMI Gal. Poles 



io- 26 



io- 26 



10 2 I0 3 
m x [GeV] 



III 4 




10 -20 rr- 



10- 21 E 



> 

o ^ - ' 



10- 



~ FERMI 3x3 
FERMI 5'x 30" 
FERMI 10"- 20" 
FERMI Gal. Poles 

10 3 















O 

= o 
= oo 










" FERMI 3"x 3' 1 

FERMI5"x30" 

FERMI 10"- 20' e 

FERMI Gal. Poles I 







m x [GeV] 
DM DM -> WW, NFW profile 











: > 




SOGi 




;i 




5 


FERMI 3"x 3" ! 

FERMI 5"x 30° 

FERMI 10'- 20" e 

FERMI Gal. Poles - - : 







io- 21 



10- 22 5 



10 10 3 
m x [GeV] 

DM DM -> WW, Iso profile 




> 
o 



FERMI 3x3 
FERMI 5'x 30' 
FERMI 10'- 20" 
FERMI Gal. Poles 



I0 2 10 3 
m x [GeV] 



10 2 10 3 
hi, [GeV] 

DM DM -> tt, NFW profile 











: > 




<u 

O 




5 O 
= OO 




;i_ 






FERMI 3" x 3' 1 




FERMI 5" x 30" 




FERMI 10"- 20' s 




FERMI Gal. Poles : 







10 2 10 3 
«i, [GeV] 

DM DM -> tt, Iso profile 



io 4 




FERMI 3x3 
FERMI 5"x30" 
FERMI 10°- 20' 
FERMI Gal. Poles 

10 2 10 3 
m x [GeV] 



IO 4 



Figure 3: Like figure 2, but for DM annihilations into bb, W + W , and tt. 



o For NFW or Einasto profiles, a significant contribution to the positron fraction appears 
strongly excluded by gamma-rays. Note that for cored-profiles, the inner-galaxy data are 
relaxed, but the higher latitude ones are not and those are sufficient to suggest an (at most) 
subleading contribution of DM to the positron flux. 

o For NFW or Einasto profiles, already this extremely conservative analysis gives constraints 
on (av) which, for typical masses m x ~ 100 GeV, are within one order of magnitude of the 
value expected for a S-wave annihilating thermal relic. 
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DM -» ee, Einasto profile 



DM -» Einasto profile 



DM -» rr, Einasto profile 




Figure 4: Similarly to figure 2 but for decaying Dark Matter. The vertical axis reports here the 
half-life Td ec in seconds. The exclusion contours are due to Fermi observations of the 10° — 20° 
strip' (red dashed line), the \b\ > 60° 'Galactic Poles' region (black long dashed line) and the 
isotropic flux (magenta dotted line). We also report the regions that allow to fit the PAMELA 
positron data (green and yellow bands, 95% and 99.999% C.L. regions) and the PAMELA positron 
+ Fermi and HESS data (red and orange blobs, 95% and 99.999% C.L. regions) in terms of decaying 
Dark Matter. We here report only the case of an Einasto galactic DM profile: the cases of an 
Isothermal or a NFW profile are essentially identical (see text for details). 



Finally, in Fig. 4 we report the exclusion plot for the case of leptonically decaying DM. We show 
that: 

o The intermediate and high latitude galactic observations impose significant constraints, that 
are in general stronger that any other Galactic gamma ray or neutrino constraint (see 
e.g. [13]). They are however insufficient to probe the PAMELA+Fermi+HESS regions. 

o The residual isotropic radiation measured by Fermi imposes the strongest constraint. It ex- 
cludes the decay explanation for the PAMELA+Fermi+HESS anomalies for the t + t~ channel 
independently on the profile. Even the allowed region for the /x 4 ' [i~ channel starts to be con- 
strained by the data and it is reasonable to expect that a more refined analysis with future 
data will be able to probe definitely this region as well. 

o Note that in Fig. 4 we report only the case of the Einasto DM profile: the exclusions plots 
for the NFW or cored isothermal cases are essentially identical. This is to be expected: the 
allowed PAMELA+Fermi+HESS regions bear a very weak dependence on the profile choice, 
as positrons and electrons come from the local halo region where all profiles resemble one 
another. The DM galactic gamma ray signal from intermediate and high latitudes (and 
therefore the constraints) vary only within 5% or less, as illustrated by the values of the Jd ec 
factors in Table 1 (the isotropic flux is obviously halo independent as already discussed). 



4 Discussion and Conclusions 

In this article, we have provided a first assessment of the power that new data on the diffuse 
emission from the Fermi Gamma-ray Space Telescope have in constraining Dark Matter indirect 
signals. Even under the very brutal approximation of neglecting any astrophysical background 
contributing to the signal and using conservatively 3 a exclusion criteria, current data from the 
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Figure 5: For one selected observational region (the 10° — 20° strips') and profile, and for two 
exemplar annihilating DM candidates, we show the current bounds derived by our conservative 
procedure of requiring that the DM signal must not exceed the observational data points by more 
than 3a (red dashed line) and the prospective regions that can be explored by Fermi if an astro- 
physical gamma ray background is assumed and the exclusion criterion is relaxed to 2a (blue short 
dashed line). 



inner Galaxy (e.g. '3° x 3°') exclude a benchmark DM mass m x ~ 100 GeV if its annihilation 
is larger than a factor 5-=-30 (depending on the channel) of the typical (av) ~ 3 x 10 _26 cm 3 /s, 
when profiles suggested by N-body simulations are employed. Higher-latitude constraints are a 
factor ~ 10 weaker and comparable to constraints for cored profiles. It is remarkable that already 
such a simplified analysis is powerful enough to explore regions of parameter space not excluded 
otherwise, providing constraints which are comparable to or better than those obtained by the 
Fermi collaboration by analyzing dwarf spheroidals [41] or the isotropic signal from cosmological 
DM annihilation [42]. This confirms, if needed, the Galactic halo as the "target of excellence" for 
constraining or detecting gamma rays from DM. 

On the other hand, the absence of astrophysical background is an extremely (unrealistically) 
conservative assumption as visual inspection of the plots in Fig. 1 confirms. In the pre-Fermi 
era, some studies have been performed showing the possible improvement in sensitivity when 
accounting for pointlike and diffuse sources in the Galactic Center region (see e.g. [43]). The 
current high-quality data certainly allow one to improve over these exploratory studies to forecast 
the ultimate Fermi sensitivity to DM. While a proper treatment of this problem goes beyond our 
current purposes, in Fig. 5 we present for illustration the exclusion plots that would follow from the 
current '10° — 20° strips' data if its bulk could be robustly attributed to astrophysical processes, 
as in the adjusted propagation model shown in [26] and the exclusion criterion is relaxed from 3a 
to 2a. The 'improvement' is about a factor of 2. Likely, intermediate-latitude DM bounds could 
be made competitive with current conservative inner-galaxy constraints. In turn, the latter could 
improve significantly if maps were cleaned from further astrophysical sources contaminating the 
total flux: notice that the '3° x 3°' degree field data are not corrected for pointlike sources [24], 
which have already been identified by the very Fermi satellite (see e.g. [44]). It appears realistic 
to expect that by the final stages of the Fermi mission a better knowledge of the astrophysical 
sources (both pointlike and diffuse) would allow for significantly stronger constraints or perhaps 
enough sensitivity for a detection of a "benchmark DM candidate" , in particular if the DM halo 
profile is steep enough toward the GC, or if significant substructures are present enhancing the 
diffuse flux. This would be essentially consistent with early expectations [45]. 
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We have also addressed the issue of the consistency of DM interpretation of the leptonic CR 
spectra from PAMELA, Fermi, HESS with the new diffuse gamma ray data. It is worth stressing 
that, even before such new Fermi data, the annihilating DM interpretation of the leptonic CR 
spectra was challenged by the constraints from gamma rays and radio waves from the galactic 
center [46, 47, 48], from the integrated cosmological flux of ICS photons [49, 50, 51], from neutrinos 
from the galactic center [52, 13] and from Big Bang Nucleosynthesis [53]. Among the most robust 
constraints are those from CMB observations and the reionization history of the Universe [54, 55, 
14, 51, 56], which rule out (if barely) the best fit regions of PAMELA+Fermi+HESS. 

The new diffuse gamma ray data allow a further test of the DM hypothesis, in particular of the 
prompt and especially IC emission at relatively large latitudes (for former studies with previous 
data see e.g. [57, 12, 13]; for related studies with current data see [58, 67]). These signals have a 
weaker dependence on the extrapolations of the profile in the inner region of the Galaxy than most 
previously cited constraints. For example, the strip 10° — 20° does not depend at all on the profile 
in the inner ~ 1.5 kpc of the Galaxy and most of the signal comes actually from a significantly 
larger distance, due to geometric effects. Moreover the new data allow for to stringently probe 
the decaying DM explanations of the anomalies. 

In summary, here are our results: 

• Decaying DM appears significantly excluded from the isotropic emission inferred from the 
Fermi data, unless it decays mostly/exclusively in fi ± . 

• We find that DM annihilating hadronically -already excluded by antiproton data as a joint 
explanation for the e + and e + + e~ anomalies- cannot even account for a significant fraction 
of the positron fraction. This remains essentially true even for cored DM profiles. 

• DM annihilating leptonically is excluded as leading explanation for CR charged lepton data 
in all cases if the halo shape is close to the Einasto or NFW one, and disfavoured even for 
cored profiles for r mode. It appears that only for cored halo profiles, mostly annihilating 
into /x 1 * 1 , some viable fit of the CR leptons can be found. However, it is reasonable that even 
such "fined-tuned" situations may run into some troubles when scrutinized more deeply. 

These conclusions are in general based on robust and conservative procedures, as discussed in the 
text. We list here a few fine points and caveats. 

In realistic models where DM annihilates into several channels -including quarks (— > 7r's) 
and gauge bosons- the B.R. into jj^ might be insufficient to evade the bounds. A similar effect 
would derive from considering three-body channels with final state electroweak radiation, which 
in particular for these engineered "leptophilic" scenarios could have non-negligible B.R. [59]. 

Also, when moving from the attitude of constraining an otherwise undetected DM to the 
hypothesis of checking predictions of a DM candidate postulated to account for some other CR 
signal, omitting astrophysical background in the gamma-channel is not a self-consistent approach. 
Indeed, in order to account for e + and e~ data, some additional background must be introduced 
and some choices for propagation parameters have been made (e.g. normalization and index of 
diffusion function, equivalent to a "grammage" in a leaky box picture). The low-energy e + fraction 
is explained by invoking CR spallation and n — > fi — > e production, etc. Thus, any "best-fit choice" 
with background plus DM is accompanied by a corresponding lower limit for the diffuse gamma- 
ray background due to CR spallation in the ISM. While in principle other sources of gamma's not 
showing-up in electrons are possible, the one above is unavoidable. It is worth reminding that 
predictions based on the "conventional" model calculated by GALPROP [60] and recently refined 
by the Fermi team on the light of their new data [26] fit the data quite satisfactorily without the 
need for a leading role of DM. A similar caveat applies to the practice of best-fitting independently 
the background for different channels: since different astrophysical backgrounds are correlated, a 
proper analysis may restrict the parameter space for DM models to fit the data. 
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We conclude this article by a few comments on how to modify/improve the constraints. 

> The HESS (e + + e~) data points might contain a sizable contamination from diffuse isotropic 
gamma rays, up to 50% [6]. Strictly speaking, for a given DM candidate it is therefore the 
sum of the (e + + e~) yield and the isotropic 7 (as computed in this paper) which has to be 
compared to the HESS data. However, for the typical candidates that we considered, the 
isotropic gamma ray flux turns out to be much smaller than the corresponding (e + + e _ ) flux. 
No more stringent bounds can therefore be derived by this procedure. For candidates with 
larger mass (above several TeV) this constraint could however turn out to be significant. 

> As already mentioned, should the Fermi collaboration present data or limits on the flux at 
the highest energies (e.g. the E ~ 300 GeV bin) they would be very powerful in constraining 
DM models that still provide marginal fits to the data, provided that the associated error 
bars prove to be small enough. Even an upper limit at the level of the current detected flux 
around ~ 100 GeV would strengthen the bounds by a factor ~2 or so. 

> Also, we have not scanned the Fermi maps looking "a posteriori" for the most constraining 
regions for a given DM model starting directly from the 7— maps. If anything, such strategy 
would strengthen the constraints. 

> Moreover, note that Fermi has presented some early results on the limits on gamma-ray 
line-emission [40]. This emission is unavoidable at loop level, but more model dependent. 
For a typical B.R.~ 10~ 3 and for relatively light candidates (up to a 200 (400) GeV for 
annihilating (decaying) DM) the resulting bounds should be comparable or stronger than 
the ones reported in the text. 

> Additional constraints are expected from the associated neutrino production (unavoidable 
already from the two-body channel for the fj, and r modes). It is likely that a reanalysis 
of the Super-K data, including large regions around the Galactic Center, may provide a 
conclusive test of these models (or, alternatively, an analysis of future Icecube+DeepCore 
data [61].) 
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A Uncertainties of the constraints due to galactic magnetic field mod- 
els and e ± synchrotron radiation 

In this appendix, we address the issue of how the constraints that we presented are affected by 
the uncertainties on the galactic magnetic field. Indeed, in our analysis we have always considered 
(in line with the analysis of [12]) the energy losses suffered by the DM-produced electron and 
positrons into synchrotron radiation as subdominant. If the magnetic field is however significantly 
large, then competing synchrotron energy losses can be important, the IC emission can be reduced 
and the propagation of e fitting the PAMELA+Fermi+HESS modified. 3 

3 We stress that of course only the IC compton emission from DM is affected by the value assumed for the 
magnetic field B. This is the dominant signal only for the leptonic modes, in particular e + e~,fi + ^i~. The prompt 
component (which is e.g. present or dominant in the t + t~ annihilation or decay channel) is not affected. 
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For E e 3> GeV and in the Thomson regime, the energy loss scales roughly as (e.g. [62], 
Eq. (4.6.1)) 

^ -E 2 [B 2 ± + 27 u 7 ] (11) 

where B± is the perpendicular component of the magnetic field with respect to the direction of 
motion, measured in fiG, and w 7 is the energy density in low-energy photons, measured in eV. In 
the vicinity of the solar system, w® « 1 in these units (see e.g. [62], Table. (2.2)), while in the GC 
region one has u GC ~ 10 (see for example Fig.l in [12]). From the relation above, we can derive 
that as long as Bf_ <C 5.2 and B GC <C 16.4, the IC losses dominate. Equivalently, for an isotropic 
field for which B± m y^2/3B to t, the IC dominates as long as 

{Bf ot < 6.4, B GG < 20.1}, (in /iG units). (12) 

On the other hand, the interest of DM annihilation/decay final states like e + e~,/i + /i~ is mostly 
related to phenomenological reasons, namely the fit of Fermi and/or PAMELA lepton data. So, 
rather than the absolute level of the bounds, it is more interesting to know how the relative position 
of the best-fit region vs. bounds moves with changing parameters. In the (physically unlikely, see 
below) limit where the magnetic loss term dominates over IC, what alters this position of, say the 
GC bounds from IC compared to the fit, are (to a first approximation) differences in the ratio 
(Bf) 2 /u® (the quantity relevant to the local propagation of the that fit PAMELA+Fermi+HESS) 
compared to (B 1 ^ ™) 2 / '■u^ here where with the notation there we want to indicate that this is intended 
evaluated in the regions associated with the window of gamma ray observation and therefore it 
is the quantity relevant to energy losses suffered by the that produce the bound. For the 
large regions in our analysis, the conditions there are very much expected to be similar to the 
conditions here, i.e. at . A larger sensitivity to the above considerations is instead obtained 
when considering the bounds coming from the 3° x 3° region around the GC. We consider already 
the bounds coming from this region as being the least robust, due to the fact the we neglect the 
effect of diffusion. 

In order to estimate the possible error introduced by the ignorance on the -B-field, first of all 
we consider the models discussed in [63]. We perform an average of B 2 over a kpc scale volume, 
whose square root we dub for simplicity of notation (B &,GG ). This is the typical propagation scale 
of high-energy electrons and roughly the size of the 3° x 3° region. For the two models (TT and 
HMR) which are regular towards the GC region, one gets values well within the requirements of 
Eq. (12), for example {(B GC ), (B G )} = {4.6, 1.4} for the TT model. 

These models however only aim at describing the large scale features of the field as inferred e.g. 
from Faraday Rotation measurements. Since additional power may be present in the smaller scale 
field to which synchrotron radiation is responding, too, a perhaps more appropriate representative 
for the total field can be taken from [66], 

B tot = B exp(-(r - R Q )/R B ) exp(-\z\/z B ) ■ (13) 

The suggested parameters i?B=10 kpc, I?o=6.1 /iG, and zb = 2 kpc reproduce the absolute 
magnitude and profiles of the 408 MHz radio emission; for these values one gets {(B GC ), (B G )} = 
{5.4, 12}, again consistent with the requirements of Eq. (12). 

Finally, one can consider models which extrapolate to a very high value towards the GC, like 
the PS one in [63], where there is a mathematical divergence at the GC position. Those kind of 
models have to be regularized based on physical arguments. In [64] it was extensively argued that 
the value of the pervasive magnetic field in the inner degrees of the Galaxy is likely ~ 10 /iG, i.e. 
well below the critical value of ~ 20 fiG reported above and consistent with the estimates from 
previous fitting formulae. On the other hand, high fields of the order ~ 100 fiG or larger might 
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be present in relatively thin filaments (see [65] for references). For the purposes of this paper, it 
is sufficient to report as a general argument that as long as 

f^Vx^«l (14) 
\2QfxGj lkpc 3 1 ; 

the limits in Eq. (12) are satisfied and the bounds derived in the main text are robust. For 
example, if filaments were to occupy a (100 pc) 3 volume, where the field reaches ~ 100 fiG field, 
the above limit is easily fulfilled. What if even the above condition is violated? Of course, the IC 
bound from the inner Galaxy would be weakened, since most of the energy loss would be damped 
into synchrotron radiation. The recent analysis in [65] however showed that, when combining 
radio and gamma-data, the overall constraints on the annihilation cross section are quite robust, 
independently of the 5-field assumed, although which channel yields the best constraint has a 
-B-dependence. 

At the light of these results, we conclude that the IC bounds are quite robust against varia- 
tions in B field, except possibly from the region very close to the GC and for extreme values of 
parameters. In the latter case, however, deriving more reliable constraints from IC is not only 
model-dependent, but actually of little relevance, since other bounds are known to be as stringent 
(or more stringent) than the IC ones. 
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